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Abstract 

We investigate the behavior of the quasi-particle energy gap near quantum phase transitions in a two- 
dimensional disordered boson Hubbard model at a commensurate filling. Via Monte Carlo simulations of 
ensembles with fixed numbers of particles, we observe the behavior of the gap as a function of the tuning 
parameter for various strength of diagonal disorder. For weak disorder, we find that gapped Mott insulating 
phase is sustained up to the transition point and disappears only in a superfluid, strongly supporting a direct 
Mott-insulator-to-superfluid transition. Bose glass behavior, insulating with vanishing gap, appears only 
when the strength of disorder is bigger than a critical value. 

PACS numbers: 74.78.-w, 73.43.Nq, 74.40.+k 
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I. INTRODUCTION 



The interplay of interaction and disorder in quantum phase transitions has been an important 
issue for past decades. In systems of bosonic particles, the phase coherence in a superfluid(SF) 
can be destroyed by either interaction or disorder effect. However, it is not obvious how these two 
effects interplay to drive the quantum phase transitions. To investigate the role of disorder, the 
superfluid-insulator transitions of the disordered boson Hubbard model 1 in low dimensions have 
been studied intensively 2 i 3 i 4 ^i 7 i 8 i 9 i 1Q i 1 12 . Experimental realization of this problem may include 
the transitions in 4 He in porous medial, cold atoms in optical lattice o 14 ' 15 , superconducting thin 
films 16 , and small Josephson-junction arrays 17 . 

In the absence of disorder, a strong interaction suppresses the density fluctuations and results in 
a Mott insulator (MI). In the phase of MI it costs finite energy to add an extra particle. This energy 
gap vanishes at the MI-SF transition as the hopping strength increases, allowing delocalized quasi- 
particles without any energy cost. In the presence of disorder, the situation is subtle. Based on 
a naive single-particle picture, it has been argued that an arbitrarily weak random potential will 
localize a single-particle excitation as soon as the Mott gap vanishes. Therefore, there exists an 
insulating phase with vanishing gap, called the Bose glass (BG) phase, which intervenes between 
MI and SF phases. A hallmark of the BG-SF transition is that the value of the dynamical critical 
exponent is given by z — d, where d is the dimensionality of the system, since the BG phase is 
compressible 1 . 

Many numerical searches of the BG phase have been performed by identifying the value of the 
critical exponents including z in two dimensions. Though the existence of BG for the strongly 
disordered case has been generally accepted, for the weakly disordered case it remains as a 
controversial problem. Direct MI-SF transitions, without the intervening BG phase, have been 
foun d 18 ' 19 ' 20 ' 21 , especially at commensurate fillings, while no clear numerical evidence of the BG 
phase has been reported in this case. This may suggest that the physics of disordered system 
is divided into the strong and the weak regimes, including the possibility that the particle-hole 
excitations may screen out weak random potentials. So far most of the numerical works focus 
on whether the critical exponents support the direct MI-SF transition. Therefore, it is very use- 
ful if one can directly observe the gap near the superfluid-insulator transition to clarify whether 
there exists the gapless insulating BG phase between the MI and the SF phases in the presence of 
disorder. 
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In this work, we investigate the behavior of the quasi-particle gap as a function of the tuning pa- 
rameter for different strength of disorder in a disordered two-dimensional boson Hubbard model. 
We find that for weak disorder the gap vanishes sharp at the transition, strongly supporting the ex- 
istence of the direct MI-SF transition. This feature disappears for strong disorder, thus supporting 
the scenario that the disordered cases are divided into the strong and the weak disorder regimes. 

H. MODEL 

Interacting bosons moving in the presence of a random potential can be described by a boson 
Hubbard model 

H bH = jJ2 n ) - E m - * E( & fe + fyi), (i) 

3 3 (ij) 

where &_,-(&!-) is the boson annihilation (creation) operator at j-th site, rij is the number operator, U 
and t are strength of the on-site interaction and the hopping amplitude between nearest neighboring 
sites (ij), respectively, and Hj denotes a diagonal random potential due to local impurities. We take 
a random potential — A <//,-< A, where A is the strength of disorder, to make its average to be 
zero for the case of a commensurate density. The phase diagram of this model with possible three 
phases, known as the MI, BG, and SF, has drawn a great deal of attention, focusing on whether the 
BG phase intervening between the MI and the SF phases exists even if the disorder is weak. 

It is believed that in the limit of large density these phases and the transition between them can 
be described by the quantum rotor model 6 

H * = ^H(-qq— ft) 2 -JH oosCft - e s ), (2) 

3 3 (ij) 

where J = tn (n is the average integer number of bosons per site), = (ij/U, and 9j is the 
phase angle of bosons at j-th site. 

For numerical simulations, it is very useful to write the partition function Z = 
J2{e} exp(— @Hqr) in terms of dual current variables so that 

vj=o 

z= e~ Hcl/K (3) 
{J} 

with the corresponding classical Hamiltonian 7 

#cl[J] = E + J L + J llr)\ ' ih^r)}, & 
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where K ~ y J/U is the tuning parameter controlling the quantum fluctuations. Here V • J = 
means the current conservation condition at each site (j, r). Note that in this representation the 
winding number in r-direction is the net number of bosons, N. In other words, N = means a 
commensurate density, and N = 1(N = — 1) represent the existence of extra quasi-particle (hole), 
etc. 

When temperature T — > 0, we expect Z ps e~^ Eo , where i?o is the ground state energy of the 
quantum Hamiltonian. Thus E is closely related with the free energy of the classical action, F c i. 
This allows us to define a quasi-particle gap 

from the free energy with a fixed number N. For finite systems of size L x L x L T , where L is 
the size in a spatial direction and L T in the temporal direction, we use L T in place of the inverse 
temperature (3 = l/T. Even though the free energy, equivalently the gap, is hard to observe 
directly, its derivative with respect to 1 /K 

- r2 ^k = r{ {Hcl)(N=1) + ^ {n= ~ 1] 

-2(tf cl >^=°)} (6) 

is easily obtainable by Monte Carlo method adopting the classical Hamiltonian in Eq. @. 

Suppose the gap vanishes at a continuous transition with the dynamical critical exponent z and 
the correlation length critical exponent v. Since we have Q g ~ \K — K c \ zu as we approach the 
transition below, we expect a diverging behavior —dVt g /dK ~ \K — K c \ zu ~ l when zv < 1, as 
schematically shown in Figure[TJ However, if the gap vanishes before the correlation diverges, the 
diverging behavior of —dVt g /dK would not appear. By observing dVL g /dK, therefore, especially 
its feature of the sharp rise and the subsequent decline, we are able to investigate the behavior of 
the gap. 

The finite-size scaling ansatz of the gap from the hyperscaling relations can be constructed as 

n g = L- z X n ((K - K c )L l '\ L T /L Z ) } (7) 
where Xn is a scaling function. By taking a derivative, we have 

- = L^ z X' n ((K - K e )L 1 '", L T /L Z ) (8) 

where X' n = dXn/d(l/K) is another scaling function. Since it is believed that the gap vanishes 
for Bose glass phase and superfluid phase, we expect dVl g /dK = in those phases. 
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FIG. 1: A schematic plot for the behavior of the energy gap, U g , and its negative derivative with respect to 
K, —dQg/dK. As gap vanishes continuously at the transition, —dQ g /dK increases sharp. 

HI. MONTE CARLO METHOD 

To investigate the behavior of the gap at a commensurate density, the expectation value of the 
energy in the equivalent classical model are calculated in three replica with N — 0, 1 and — 1 
for the same complexions of the random potential fa, and, then, average over different random 
potentials are taken. We adopt a worm algorithm 22 to update the current configurations in Eq. @ 
while rejecting any updates changing the the winding numbers in r-direction in order to keep the 
number of bosons fixed. In addition, to create a pair of particle and hole, which also does not 
change the net winding number in r-direction, a conventional method of the global loop-current 
update 7 is performed in r-direction. 

Figure |2| shows the behavior of the gap as a function of K for pure case (A = 0). Near the 
critical point K c = 0.333 22 marked by an arrow in the figure, a stiff rise and the subsequent decline 
of dVLg/ d(l/ K) appears as shown in Fig. 2a. This feature becomes sharper in bigger systems, 
indicating that the gap indeed vanishes at the critical point. Interestingly, this behavior satisfies the 
scaling ansatz of the gap, given in Eq. ® as we can find in Fig. 2b. Here we just confirm that the 
known values of the exponent z — 1.0 and v = 0.67 yield a high quality scaling behavior, such as 
the crossing behavior of L z - l l u dtt g /d{l/K) at K = K c for different sizes as a function of K and 
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FIG. 2: The behavior of the gap as a function of K for the pure case, (a) A sharp decline of —dQ g /dK is 
observed near the critical point K c = 0.333, marked by an arrow, which is a signal of the vanishing gap. (b) 
The properties of the gap shows a finite-size scaling behavior with the dynamical critical exponent z = 1.0 
and the correlation length critical exponent u = 0.67. Inset shows this scaling behavior as a function of 
the scaling variable (K — K c )L l / u . High quality data collapsing onto a curve for different sizes strongly 
suggests that the gap indeed follows a scaling relation. 

data collapsing as a function of (K — K^L 1 ^ . 

The success in understanding the critical properties for the pure case through the behavior of 
the gap motivates a extension of the study to the disordered case. Figure |3] shows the behavior of 
the gap as a function of K for different A. The feature of the sharp rise and decline in dVt g /d{l/ K) 
near the critical point appears for small A. However, this feature is absent for A > A c « 0.4. 
In this strong disordered regime, we have dVt g /dK fa near the transition point, suggesting 
that only in this strong disorder regime BG phase intervenes between the MI and the SF phases. 
These results are well consistent with the previous studie o 20 ! 21 of this model through the finite-size 
scaling properties of the superfluid stiffness and the correlation functions, showing the different 
critical exponents in the strongly and weakly disordered regimes separated by the critical strength 
of disorder A c rs 0.4. 

Figure H] shows the size dependence of the gap feature near A c « 0.4. For A < A c , the feature 
of the sharp rise becomes apparent in bigger systems. No such feature appears for A > A c . 
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FIG. 3: The behavior of the gap as a function of K for different A in the systems of the size L x L x L T = 
24 x 24 x 24. The critical points for given value of A, estimated from the finite-size scaling behavior of the 
superfluid stiffness, are marked by the arrows. 

This implies that there is a qualitative change of the gap feature across the critical strength of the 
disorder, A c . 

IV. SUMMARY 

In summary, we have investigated the behavior of the quasi-particle gap as a function of the 
tuning parameter K for different strength of disorder A in a two-dimensional boson Hubbard 
model, focusing on whether the BG phase exists between the MI and the SF phases in the presence 
of random potentials. The gap properties are studied in an equivalent classical model with a 
fixed number of particles. In the pure case, the properties of the gap satisfies a finite-size scaling 
relation with a sharp signal indicating that the gap vanishes at the transition. In the disordered 
case, our results strongly support a scenario that the cases are divided into the strong and the 
weak disordered regimes. For weak disorder, we can find the features of the vanishing gap at the 
transition, implying that the direct MI-SF transition occurs. This feature does not appear in the 
strong disorder regime, which supports the existence of the compressible BG phase in this case. 
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FIG. 4: The behavior of the gap for different sizes near A c . A qualitative change of the gap behavior 
appears across A c « 0.4. 
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